
selfefficacy<-c(0.87, 0.96)
selfefficacy.se<-c(0.03, 0.017)
selfsacrifice<-c(0.87,0.88)
selfsacrifice.se<-c(0.03, 0.026)


var.names<-c("Baseline", "Endline")
x.axis<-length(var.names):1
adjust<-.05

par(mar=c(5,5,2,2),lheight=.5)
plot(x.axis-adjust, selfefficacy, type="p", axes=F, xlab="",ylab="Share Willing to Participate in Protest",pch=19,col="black",cex=0.75,ylim=c(0.5,1),xlim=c(max(x.axis+.1), min(x.axis-.1)),main="")
axis(2,pretty(c(0,1),10), cex.axis=1)
axis(1,at=x.axis,label=var.names,las=1,tick=T,cex.axis=1)
box=(bty="1")
segments(x.axis-adjust,selfefficacy-qnorm(0.95)*selfefficacy.se, x.axis-adjust, selfefficacy+qnorm(.95)*selfefficacy.se, lwd=1.3)

segments(x.axis,selfsacrifice-qnorm(.95)*selfsacrifice.se, x.axis, selfsacrifice+qnorm(.95)*selfsacrifice.se, lwd=1.2)
points(x.axis,selfsacrifice, pch=21,cex=1,bg="white")


#add legend
legend(2,.7, c("Power of Faith Treatment Group", "Obligation Treatment Group"), pch=c(19,21),col=c("black","black"), bty="n", cex=1)

